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We propose a quantum oscillation experiment by which the rotation of an underdoped 
YBa2Cu306+a; sample about two different axes with respect to the orientation of the magnetic 
field can be used to infer the shape of the in-plane cross-section of corrugated Fermi surface cylin- 
der^). Deep corrugations in the Fermi surface are expected to give rise to nodes in the quantum 
oscillation amplitude that depend on the magnitude and orientation of the magnetic induction B. 
Because the symmetry of electron and hole cyclinders within the Brillouin zone are expected to be 
very different, the topology can provide essential clues as to the broken symmetry responsible for 
the observed oscillations. 



The recent discovery of magnetic quantum oscillations 
in high temperature superconductors provides an entirely 
new perspective on their electronic structure [TJ [5J [3J |H 
[5]. These experiments provide evidence for a Fermi sur- 
face and possible evidence for a broken symmetry ground- 
state in underdoped samples, yet the precise form of the 
symmetry remains the subject of debate. A key point 
of discussion is the carrier sign of the pocket respon- 
sible for the largest amplitude quantum oscillations in 
YBa2Cu306+a; and YBa2Cu40s [TJ E]. Quantum oscil- 
lation experiments cannot directly infer the sign of the 
carriers, leaving open many possibilities for the interpre- 
tation of the negative Hall coefficient and topology of the 
Fermi surface El El El OH EE2 EH • The choice between 
electron and hole pockets cannot be taken lightly. If the 
carriers are truly electron-like then this places limitations 
on the types of theoretical model that can be applied to 
the underdoped metallic state of high T c superconduc- 
tors. Given the possible implications for the competing 
order parameters in YBa2Cu306+a;, a thorough experi- 
mental investigation of the Fermi surface topology is of 
paramount importance. 

In this paper we show that crucial clues as to the nature 
of the broken symmetry phase that competes with super- 
conductivity may be found by performing quantum oscil- 
lation experiments in which the crystalline c axis of the 
sample is tilted away from axis of the magnetic induction 
B in different directions. The recently proposed deeply 
corrugated form of the Fermi surface cylinder (s) [HI [15] 
implies the existence of beats in the quantum oscilla- 
tion amplitude as a function of the magnitude and ori- 
entation of B. By tracking the position of nodes in the 
quantum oscillation amplitude, the shape of the in-plane 
cross-section of the cylinders can be mapped out in a sim- 
ilar fashion to angle-dependent magnetoresistance oscil- 
lations (AMRO) [TB] • A distinction between electron and 
hole pockets can be made because they are predicted to 
have different symmetry properties within the Brillouin 
zone (as depicted for the simple case of Q = (tt, tt) order- 
ing in Fig. [T^i), the details of which may further depend 
on the precise form of the broken symmetry. 




FIG. 1: a Schematic Fermi surface for the simple case of 
Q = (tt, tt) density-wave ordering, leading to reconstructed 
electron (pink) and highly elliptical hole (cyan) pockets. The 
electron pockets are expected to vanish in strong-coupling 
models, b Schematic warped Fermi surface with an elliptical 
cross-section with its major axis at an arbitrary orientation 
£ within the k x , k y plane. The extremal cross sections or- 
thogonal to k z and c are depicted in green. The red curve 
shows an example orbit when the magnetic induction B is 
rotated by a 'zenithal' angle 9 away from the k z axis. The 
half caliper width of the in-plane cross-section k\\ depends on 
the 'azimuthal' component <f> of the magnetic field orientation 
and on £. 

Independent of the broken symmetry responsible for 
the creation of a given corrugated Fermi surface cylinder, 
the &; z -dependence of the cross-sectional area Af. can be 
expanded as a series of harmonics. A schematic of a 
corrugated Fermi surface cylinder is shown in Fig. [T]d. If 
the in-plane crystal lattice vectors a and b are orthogonal 
to the interlayer vector c, A^ can be considered as the 
sum of even harmonics 

A k = A n cos(nk z c) (1) 

n 

where c = |c|. Here, Aq is the mean cross-sectional area 
of the Fermi surface, A n >i ss 47rm*i Ci „/fi 2 represent har- 
monics of the corrugation and m* is the in-plane effective 
mass. In the presence of a magnetic induction B, the k- 
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space area Ag t< p of a given semiclassical orbit (depicted in 
red in Fig.[T]D) depends on the 'zenithal' 8 and 'azimuthal' 
4> orientation of B with respect to the corrugated cylin- 
der (as depicted in Fig. [TJj). The orbit area in such a 
situation was found by Yamaji |17j to have the form 

Ae^cosO = ^yl„ J (nfcyc tan 6>)cos(nfc Zi0 c)- (2) 



Its further dependence on the half caliper width fcii of 
the corrugated cylinder within the k x ,k y plane was con- 
sidered in subsequent works [16 . Here Jo(^) is the zeroth 
order Bessel function, while fc z ,o is the k z coordinate of 
the geometric center of the orbit. If the in- plane cross- 
section is elliptical in form, as if often approximately the 
case for small pockets of carriers [16] , then 



fen = fcpJ ?ycos 



>-C) + -sin 2 (0-C), (3) 
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where k-p = JA~opK is the effective Fermi radius, rj is the 
cllipticity of the pocket (or the ratio of the minor axis 
to the semi-latus rectum, or the ratio bsetween major 
and minor axes) and £ is the orientation of its major axis 
with respect to k x within the k x , k y plane. For rj = 1, this 
reduces to a perfect circle where hi = fcp, as was origi- 
nally considered by Yamaji [IT]. The k z ^-dependence of 
the quantum oscillation frequency is obtained applying 
Onsager's relation 
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to Eqn ([2]) . Finally, the form of the magnetic quantum 
oscillations is obtained by integrating over k z ,o- On ne- 
glecting oscillations of the chemical potential (assuming 
that the corrugation is large), the oscillatory magnetiza- 
tion is given by |18j 
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are the conventional amplitude-modifying prefactors and 
X. p = 1 / K 2 "pra*k\\T jheB 1 18] . Here, p is the quantum 
oscillation harmonic index. The term containing the hy- 
perbolic sine is the thermal damping factor, the expo- 
nential factor accounts for quasiparticle scattering, while 
the cosine term accounts for Zeeman splitting (or 'spin 
splitting') [18]. 

To understand the basic form of the B-dependcnt 
quantum oscillation amplitude, it is instructive to first 
consider the simpler situation in which the second har- 
monic of the corrugation A n >2 is neglected (along with 



higher harmonics). In this case, the integral given by 
Eqn ([5]) has the simple solution 

f2irpAF e ,i 
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where AFg^ = AFJo(fc)|ctan6')/cos6' is the difference in 
frequency between maxima ('belly') and minima ('neck') 
frequencies of the corrugated cylinder (with AF being 
that at 8 — 0). In AMRO experiments, one is typically 
concerned with the 'magic' angles at which AFg^ van- 
ishes |161 117] — where the electronic structure mimics that 
of an ideal two-dimensional metal. However, nodes in the 
quantum oscillation amplitude, which occur for a finite 
value of AFg ^, can be determined with greater accu- 
racy in quantum oscillation experiments. Since AFg^ > 
B/8tt in the limit where we consider nodes, Equation ^ 
can be accurately approximated by the superposition 
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of neck F DCck = F - AFg i<f> and belly F bc n y = F + 
AFg^ quantum oscillations, giving rise to a distinctive 
beat pattern. Nodes occur whenever the -Fneck and -Fb e ii y 
oscillations are ir out of phase, which occurs at a different 
set of 'magic' angles 8 m given by 



Jo(/c||ctan0 m ) = 



(4m + 3)B 
8AF 



cos &„ 



(9) 



where m is an integer. The monotonic dependence of 
the m th node 8 m on k^/kp in Fig. 2 (in which m = 1 
and are plotted) implies that its <7>aependent angular 
position can be used to map the shape of the cylinder 
cross-section. 

We now turn to ^-dependent simulations of M z in un- 
derdoped YBa2Cu306+a; for several different Fermi sur- 
face topology scenarios. Recent quantum oscillation mea- 
surements on YBa 2 Cu3 6 . 5 4 [TH [TS] (detected over a 
wide interval in magnetic field B w P-qH) provide es- 
timates for the relevant parameters. With the exception 
of the weak 637 T frequency [2], the full waveform from 
Ref. |14j can be reproduced in Fig. [3] by numerical in- 
tegration of Eqn. Q and the adjustment of only four 
parameters [TI5] (listed in the figure caption). 

We begin by considering scenarios in which the 
F ~ 500 T frequency corresponds to an electron pocket 
cylinder, for which a single value of fcy is expected at each 
azimuthal angle cf>. Several different proposals for the 
groundstate of YBa2Cu 3 6+I predict electron pockets 
that are approximately circular in cross-section and situ- 
ated near k = (n, 0) in Fig. [TJl. These include Q = (tt, tt) 
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FIG. 2: A plot of cot6* m versus fc||/fcF according to Eqn Q 
for m = and m = 1, using c = 11.68 A and AF/B = 1 
(approximately applicable to YBa2Cu30e.54 at B — 40 T). 
For simplicity, we consider only the nodes occurring when 
fc|ictan0 m < & 
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FIG. 3: Simulation of the magnetization (or magnetic torque 
t oc —M Z B) from Ref. [TJ] using Eqn. Q and four ad- 
justable parameters Fo = (h/2ire)Ao ~ 513 T, AF = 
(h/2ne)A 1 w 40 T and AF 2 = (fi/27re)A 2 « -10 T in Eqn. Q 
and T = 170 T in Eqn. These values yield F llcck « 463 T 
and Fbdiy ~ 543 T. We have assumed the phase to be a fixed 
quantity in Eqn. (J5| [TS]. We have also assumed T = 1 K, 
m*/m c = 1.7 and (? = to be fixed quantities, although inac- 
curacies in T or m* are easily compensated by an adjustment 
of r and thus have little affect on the waveform. 



spin- or d-density wave ordering (4[ [10] , certain types of 
spin stripe ordering in which Q = (7r[l ± 25], it) where 
<5 = | and helical or spiral ordering jU [TH] . Strong 
coupling treatments in which the hopping of quasipar- 
ticles is frustrated by antiferromagetic correlations can 



also predict circular electron pockets under certain con- 
ditions [11], although in this case they are situated at 
k = (it, it). In all these cases, a circular electron pocket 
implies that fcy sa k-p for all </>, in which case nodes corre- 
sponding to m = 1 and should be observed at Q\ w 40° 
and 6*o » 55° when B = 40 T at all values of </> (as shown 
in Fig. [iji). Figure [5^ shows the calculated ^-dependent 
magnetization (using the same parameters as in Fig. |3| . 
Nodes at 81 and #0 are clearly evident in the simulation 
made using Eqn. (JsJ) in which the second harmonic of the 
corrugation A2 is neglected (black curve), and remain dis- 
cernible and at the same positions in the simulation made 
using Eqn. ([s]) in which A2 is included (green curve). 

Under certain conditions, the electron cross-sectional 
area can deviate from an ideal circular form. In S = | 
stripe models that incorporate both spin and charge or- 
dering (e.g. Figs. 3c and 3d from Ref. [7]), for example, 
the cross-section is expected to become elongated along 
k x . In Fig. [4]) this situation is approximated by an el- 
lipsoidal pocket with rj — 2 and £ = 0. The reduced 
symmetry of the pocket is reflected in the ^dependence 
of the angular sweeps in Fig. [5^, and b. For <p — 90° (in 
which, since C = 0, the axis of rotation corresponds to 
the major axis) Q\ and #0 are pushed out to very high 
angles in Fig. [5Jd, becoming practically unobservable for 
the simulation made using Eqn. ([5]) in which A 2 is in- 
cluded. For <j> = 0, by contrast, 9\ and 9 are reduced 
to small angles in Figs|4]3 and [5};. If the unreconstructed 
Fermi surface is more circular in form |10j . spiral or d- 
density wave order can give rise to trapezoidal-shaped 
electron pocket cross-sections, as shown in Fig. |4j: (using 
Q = (tt[1 + 25], it) where 5 = 0.1. In spite of the loss 
of rotational symmetry, the ^-dependent fcy is similar to 
that of an ellipse, but with its major axis along k y . Since 
the (f> dependences of 0\ and 82 are rotated by 90° with 
respect to Fig. |4]d, their measurement can by used to dis- 
tinguish stripe and spiral or d-density wave scenarios, or 
perhaps other scenarios [H] 

Turning now to the case where the F ~ 500 T fre- 
quency originates from hole pockets [TJ [HI d2], the in- 
dependences of fc||, 6\ and #0 m Fig- pi and the forms 
of the ^-dependent curves in Figs. [5}i-fbecome very dif- 
ferent from those of electron pockets. In the hole pocket 
simulations we have assumed 77 = 3.5, which is approx- 
imately what might be expected were the Fermi surface 
obtained in recent angle-resolved photoemission experi- 
ments [9] folded about the (it, 0) — (0, it) line. Since there 
are two different orientations ( sa ± 45° to be counted 
in Fig. EJi and Eqn. (|5j), two different values of k\\ give 
rise to four simultaneously beating frequencies, leading to 
more complicated ^-dependent curves in Figs. |5}i-f. On 
the other hand, the ^-dependent curves exhibit a distinc- 
tive fourfold symmetry that changes rapidly with (f> owing 
to the highly elliptical cross-sections of the pockets. For 
most arbitrary values of <fi (e.g. Fig. [5J3) , the contribu- 
tion from two very different values of h\ make the nodes 



4 




FIG. 4: Polar plots (cyan) of ku versus <f> (left-hand-side) for 
different cylinder cross-sections (black), together with polar 
plots of 81 (blue) and 9o (red) versus <j> (right-hand-side). The 
estimates of 9i and #0 are made at B — 40 T, corresponding to 
AF/B = 1 in Eqn. Q and Fig. [2] Polar plot a corresponds to 
a simple circular corrugated cylinder cross-section for which 
fcii = fcp, b corresponds an elliptical electron pocket with an 

ellipticity of T] = 2 and orientation £ = predicted in | th 
stripe models that combine spin and charge ordering, while 
c corresponds to an electron pocket of unusual cross-section 
predicted by certain d-density wave or spiral ordering models, 
d corresponds to highly elliptic hole pockets (77 = 3.5) for 
which two different orientations £ « ± 45° within the k x , k y 
plane are expected, giving rise to two different values of ku 
for each value of 6. 



somewhat difficult to locate. Distinct nodes do appear, 
however, at <f> — and 90° in Fig.[5}i when the two values 
of fcii become degenerate, or near (j> w ± 15° and ± 75° 
in Fig. [5f due to a serendipitous degeneracy of d\ and 
#o- For this reason, rather detailed ^-dependent mea- 
surements are required to pin down the precise geometry 
of the cross-section of hole pockets. 

In summary, we have shown that nodes in the quan- 
tum oscillation amplitude detected in dual axis rota- 
tion experiments can be used to infer the cross-sectional 
shape of corrugated Fermi surface cylinders in under- 
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FIG. 5: Examples of ^-dependent simulations at B — 40 T 
made using Equation (|5| using the parameters listed in Fig. [3] 
caption (green curves) and using Equation |7|) in which A2 
is neglected (black curves) for several different scenarios de- 
picted in Fig. [4] Curve a is the expected ^-dependence for a 
simple circular cross-section, while b and c are the expected 
curves obtained on rotating about the major and minor axes 
(respectively) of the elliptical pocket with r\ = 2 shown in 
Fig.[4jD. d, e and f are examples of the ^-dependence expected 
for the highly elliptical hole pockets in Fig. |4ji in which there 
are two different values of ku for each value of <f> (with the ex- 
ception of <j> = or 90° (d) where they become degenerate). 
The red and blue vertical lines point to the locations of #0 
and 81 respectively. 



doped YBa2Cu306+x, which can help to answer the ques- 
tion as to whether the prominent ~ 500 T frequency orig- 
inates from electron- or hole-like carriers, as well as pro- 
viding essential clues for the groundstate description. An 
advantage of this proposed method over AMRO exper- 
iments [16j [17] is that the shape of the cross-section is 
linked directly to the pocket responsible for the quan- 
tum oscillations (which needs always to be assumed in 
AMRO experimenta) . A further advantage is that the 
predicted nodes occur at small values of 9 where H C 2 is 
lower, thereby proving less of an obstacle for Fermi sur- 
face measurements. 

There are, however, potential caveats. In the present 
simulations we have assumed that the corrugation of the 
Fermi cylinders is caused by simple interlayer hopping 



terms t c<n . More complicated forms of warping that vary 
with k x and k y may lead to some degree of variation from 
the present simulations. A A: z -dependent damping factor, 
meanwhile (for which there is presently no evidence |19l ) , 
will diminish the presence of nodes, while Zccman split- 
ting effects (presently neglected by setting g = [4 ) 
will lead to additional nodes unrelated to the corruga- 
tion. Nodes due to Zeeman splitting (often called 'spin 
splitting zeroes') typically depend only on 9, in contrast 
to those caused by corrugation described here, which de- 
pend on B, <f> and 9. 
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